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We present a reanalysis of nonstandard neutrino-down-quark interactions of electron and tau 
neutrinos using solar, reactor and accelerator data. In addition updating the analysis by including 
new solar data from SNO phase III and Borexino, as well as new KamLAND data and solar fluxes, 
a key role is played in our analysis by the combination of these results with the CHARM data. The 
latter allows us to better constrain the axial and axial-vector electron and tau-neutrino nonstandard 
interaction parameters characterizing the deviations from the Standard Model predictions. 
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PACS numbers: 13.15.+g,14.60.St,12.20.Fv 
I. INTRODUCTION 

Current solar neutrino data [1-14], in conjunction with 
reactor data from the KamLAND experiment [15] shows 
that the neutrino oscillation mechanism is the correct pic- 
ture to explain the solar neutrino physics. Solar neutrino 
experiments are also sensitive to matter effects [16, 17], 
and the combination of both solar and KamLAND data 
determines the so called Large Mixing Angle (LMA) solu- 
tion as the correct explanation to the data. For example, 
the LMA solution is quite robust against possible un- 
certainties in solar physics, such as noise density fluctua- 
tions originated by radiative zone magnetic fields [18-25]. 
Likewise, the LMA solution is also stable with respect 
to the possible existence of sizeable convective zone mag- 
netic fields [26, 27], that could induce spin-flavor neutrino 
conversions [28, 29]. In all these cases, the KamLAND 
data play a crucial role in establishing that nonstandard 
effects can only play a subleading role [30], their ampli- 
tude being effectively constrained. 

However, while constrained by the solar and Kam- 



*Electronic address: franesfeOalumni.uv.es 
^Electronic address: Omar .MirandaOf is . cinvestav.mx 
* Electronic address: mariam.tortolaOdesy.de 

^Electronic address: valleOif ic .uv. es ,URL:http: //ahep .uv . es/ 



LAND data in an important way, neutrino nonstandard 
interactions (NSI) still provide an important exception 
to the robustness of the neutrino oscillation interpreta- 
tion [31, 32]. Indeed, it has been found that they might 
even shift the solution to the so-called dark side region 
of the neutrino parameter space [33] . 

Given the envisaged precision expected in upcoming 
oscillation studies [34] , one needs to further scrutinize the 
possible role of NSI [35] . The intrinsic importance of NSI 
stems from the fact that they are characteristic features 
of theories of neutrino mass [36] and that their magni- 
tude provides important guidance in order to distinguish 
the simplest high-scale seesaw models [37-41] from those 
seesaw scenarios based at low-scale physics, such as the 
inverse [42, 43] or linear seesaw mechanisms [44], as well 
as radiative models of neutrino mass [45-47]. 

In this paper we reanalyze the robustness of the os- 
cillation interpretation of the solar neutrino data in the 
presence of nonstandard interactions. Besides all the so- 
lar neutrino data used in our previous study [33] , here we 
take into account new solar data from SNO phase III [11], 
the first real-time measurements of the solar Beryllium 
flux at Borexino [12], as well as the new and more pre- 
cise KamLAND data [15]. We have also considered in our 
calculations the new solar fluxes and uncertainties from 
the updated Standard Solar Model (SSM) [48]. We show 
explicitly that the degenerate solution in the dark side 
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region still remains plausible even after inclusion of these 
new data. Besides updating the analysis of nonstandard 
neutrino-down-quark interactions, we stress the key role 
is played by the combination of these results with the 
measurements of the electron-neutrino-quark cross sec- 
tion at the CHARM accelerator experiment. Although 
it is sensitive only to the interactions of electron neu- 
trinos, when combined with solar and KamLAND data, 
the latter allows us to improve the determination of the 
tau-neutrino nonstandard axial and vector couplings 

In what follows we will focus on nonstandard interac- 
tions that can be parametrized with the effective low- 
energy neutral currents four-fermion operator 1 : 

Cnsi = ~e f a P p 2V2G F (Va^Lvp) {f'fPf) , (1) 

where P = L, R and / is a first generation fermion: 
e,U,d. The coefficients ej^ denote the strength of the 
NSI between the neutrinos of flavours a and (3 and the 
P-handed component of the fermion /. For definiteness, 
we take for / the down-type quark. However, one can also 
consider the presence of NSI with electrons and up and 
down quarks simultaneously. Current and expected lim- 
its for the case of NSI with electrons have been reported 
in the literature [50-52]. Here we confine ourselves to 
NSI couplings involving only electron and tau neutrinos. 
This approximation is in principle justified in view of the 
somewhat stronger constraints on interactions, for a 
discussion see Ref. [53-55]. 

Nonstandard interactions may in principle affect neu- 
trino propagation properties in matter as well as detec- 
tion cross sections and in certain cases they can also mod- 
ify the assumed initial fluxes 2 . NSI effects in neutrino 
propagation affect the analysis of data from solar neu- 
trino experiments and to some extent also KamLAND, 
through the vectorial NSI couplings e^i = On 
the other hand detection shows sensitivity also to the ax- 
ial NSI couplings = s^—e^j in the SNO experiment. 
These points will be analyzed in detail in Section II, af- 
ter a brief discussion of the experimental data included 
in our study. In Section III we will focus on the study 



1 A recent study of CC non-standard interactions has been given 
in Ref. [49]. 

2 We assume a class of models of neutrino mass where NSI leave 
the solar and reactor neutrino fluxes unaffected. 



of the non-universal nonstandard interactions, combining 
the results of the CHARM experiment together with our 
solar analysis in order to obtain a new constraint for the 
tau neutrino nonstandard-intcraction with d-type quark. 
Finally we will conclude in Section IV. 

II. SENSITIVITY OF SOLAR AND KAMLAND 
DATA TO NSI 

Here we will adopt the simplest approximate two- 
neutrino picture, which is justified in view of the strin- 
gent limit on #13 [56] that follows mainly from reactor 
neutrino experiments [57]. 

A. Data 

In this subsection we will describe the input data re- 
quired to analyze the sensitivity of solar and KamLAND 
neutrino data to the presence of NSI. This will include 
not only the experimental data samples by all the de- 
tectors considered, but also the theoretical predictions 
required to simulate the solar neutrino production pre- 
scribed by the SSM. 

Concerning the solar neutrino data, we have included 
in our analysis the most recent results from the ra- 
diochemical experiments Homestake [I], SAGE [2] and 
GALLEX/GNO [3, 4] , the zenith-spectra data set from 
Super-Kamiokande I [5, 6], as well as the results from 
the two first phases of the SNO experiment [7-10]. The 
main updates with respect to our previous work [33] is the 
inclusion of the data from the third phase of the SNO ex- 
periment [11], where 3 He proportional counters have been 
used to measure the neutral current (NC) component of 
the solar neutrino flux, and the latest measurement of 
the 7 Be solar neutrino rate performed by the Borexino 
collaboration [12, 13]. 

In our analysis we use the solar neutrino fluxes and 
uncertainties given by the latest version of the SSM [48] . 
The latter provides an improved determination of the 
neutrino flux uncertainties, mainly thanks to the im- 
proved accuracy on the 3 Hc- 4 Hc cross section measure- 
ment and to the reduced systematic uncertainties in the 
determination of the surface composition of the Sun. In 
Ref. [48] two different solar model calculations are pre- 
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sented, corresponding to two different measurements of 
the solar metal abundances. For our analysis we have 
chosen the model corresponding to a higher solar metal- 
licity, BPS08(GS), although we have checked that the 
use of the lower metallicity model BPS08(AGS) does not 
change our results substantially. 

The KamLAND experiment observes the disappear- 
ance of reactor antineutrinos over an average distance of 
180 km. Given that, in their way to the detector, reactor 
neutrinos can only traverse the most superficial layers of 
the Earth, the resulting Earth matter effects are almost 
negligible. The same applies also to the nonstandard in- 
teractions we are considering. However, for consistency 
with our analysis of solar neutrino data, in our simula- 
tion of the KamLAND experiment, we have included the 
effect of NSI over the antineutrino propagation. In par- 
ticular, we have considered that neutrinos travel through 
a layer of constant matter density equal to the terrestrial 
crust density (~ 2.6 g-cm~ 3 ). Here we have used the 
latest data release from the KamLAND reactor experi- 
ment [15], with a total exposure of 2881 ton-yr, which 
brings in a big statistical improvement with respect to 
the previous data reported by the Collaboration [58] . We 
have restricted our analysis to the energy range above 
2.6 MeV where the contributions from geo-neutrinos is 
less important. As we will see in the next section, the 
inclusion of the new KamLAND data will be very impor- 
tant for the improvement of the results, given the good 
precision achieved in the determination of the oscillation 
neutrino parameters. 

B. Effects in neutrino propagation 

We first reanalyse the determination of the oscillation 
parameters in the presence of nonstandard interactions. 
The Hamiltonian describing solar neutrino evolution in 
the presence of NSI contains, in addition to the standard 
oscillations term, 
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a term -Hnsi , accounting for an effective potential induced 
by the NSI with matter, which may be written as: 



H ml = V2G F N d 
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e e 



(3) 



Here e and e' are two effective parameters that, accord- 
ing to the current bounds discussed above (e££ ~ 0), are 
related with the vectorial couplings which affect the neu- 
trino propagation by 3 : 



dV 
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The quantity N d in Eq. (3) is the number density of the 
down-type quark along the neutrino path, and 023 is the 
atmospheric neutrino mixing angle. 

From Eqs. (2) and (3) one sees that the solar neutrino 
mixing angle in the presence of nonstandard interactions 
is given by the following expression: 

Am 2 cos 26 - 2^2 EG F {N e - e 1 N d ) ... 
cos20„, = — — - , (5) 



[Am 2 



where 



[Am 2 ] 2 matter = [Am 2 cos 26 - 2^2 EG F (N e - e'N d ) 



Am 2 sin 29 + ^V2eEG F N d 



(G) 



Therefore, and as discussed in Ref. [33], there exists a 
degeneracy between the non-universal coupling e' and the 
neutrino mixing angle 6 which makes possible to explain 
the solar neutrino data for values of the vacuum mixing 
angle in the dark side (6 > 7r/4), for large enough values 
of e': 



e > 



2^/2EG F N e + Am 2 | cos20| 
2V2 EG F N d 



(7) 



For instance, for the typical values of the solar neutrino 
energies and matter densities one has e' > 0.6. Indeed, 
as we showed in [33] , the effect of NSI on solar neutrino 
propagation implies the presence of an additional LMA- 
D solution whose status we now reanalyse in the light of 
new data. 

We now turn to the combined solar + KamLAND 
analysis. Following the considerations above, we have 
performed a new analysis of all the solar neutrino data 
discussed in Section II A combined with the recent Kam- 
LAND result [15]. The main result is shown in Fig. 1. 
There, we plot the allowed regions (90, 95 and 99 % 
C.L.) in the solar neutrino oscillation parameter space 



3 For the derivation of the effective couplings in the general three- 
neutrino framework sec Ref. [59] . 




Figure 1: The left-bottom panel indicates 90%, 95% and 99% C.L. allowed region from the solar and solar + KamLAND 
combined analysis, the result of which is presented as a zoom in the top-right panel, where the shaded regions show the result 
of the current analysis while the lines indicate the earlier regions [33] for comparison. One sees that the dark side solution is 
still allowed in the presence of nonstandard interactions. The other two panels indicate the corresponding \ 2 projections. In 
all cases, we marginalize over the NSI parameters s and e . 



(sin 2 #sol, ^ m soh) obtained in the analysis of solar and 
solar + KamLAND neutrino data, after marginalizing 
over the NSI parameters in our 4-dimensional x 2 anal- 
ysis: x 2 ( sm2 #sol, A m soL> e, e'). The Ax 2 profiles as 
a function of each parameter are also shown. One can 
see that the region in the so-called dark side of the neu- 
trino parameters [33] remains even after the inclusion of 
the new data. Note however that its status is somewhat 
worse than previously. In contrast, as seen in the fig- 
ure, the other solutions LMA-0 and LMA-II [33] which 



were present before have disappeared as a result of the 
new KamLAND data, that now provide a very precise 
measurement of Am 2 . 

In order to better understand the results obtained, we 
plot at Fig. 2 the neutrino survival probabilities for differ- 
ent reference points. First we have considered the global 
best fit point from the combined solar + KamLAND anal- 
ysis, labeled as "global best fit" in the figure, with the 
following parameter values: (sin 2 #sol, Am| OL , e, e') = 
(0.32, 7.9 x 10" 5 cV 2 , -0.15, -0.10). We have also consid- 
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Figure 2: Neutrino survival probabilities averaged over the 8 B neutrino production region for three reference points (see text 
for details). Lines are compared with the experimental rates for the pp neutrino flux (from the combination of all solar 
experiments), the 0.862 MeV 7 Be line (from Borexino) and two estimated values of the 8 B neutrino flux from Borexino and 
SNO (third phase). Here, vertical error bars correspond to the experimental errors, while the horizontal ones indicate the 
energy range observed at the experiment. 



ered the best fit point in the absence of NSI: (0.30, 7.9 x 
10" 5 eV 2 , 0.00 ,0.00), labeled as "BF without NSI", and 
allowed with a A\ 2 = 2.7, and finally the best fit point in 
the "dark side" of the oscillation parameters (labeled as 
"BF dark side") with (0.70, 7.9x 10~ 5 eV 2 , -0.15, 0.95) 
and Ax 2 = 2.9. As we see all points are in perfect agree- 
ment with the low energy (pp and 7 Be) measurements. 
Concerning the most energetic boron neutrinos, where 
matter effects are more important, and as a result there 
is a higher NSI-scnsitivity of the corresponding profiles, 
the presence of NSI provides an slightly better agreement 
with the data than the standard one without NSI, mainly 
thanks to the flatter spectrum predicted above 5 MeV. 
The "dark-side" solution also gives predictions for the 
survival probability of boron neutrinos which are com- 
patible with the experimental results. From the different 
predictions obtained for these 3 reference points above a 
few MeV, one sees that this region could be crucial in or- 
der to break the degeneracy among the various solutions. 
Therefore, a better measurement of the boron neutrino 
flux with a lower threshold (like the ones expected from 
Super-K-III and SNO [60, 61]) will be of great help. On 
the other hand, a very precise measurement of the pep 
and beryllium neutrino fluxes may also contribute to lift 
the degeneracy between the standard and dark-side so- 
lutions. 



By analysing the goodness of the neutrino oscillation 
solutions in the presence of NSI one can also constrain 
the NSI parameters e and s'. In order to do this we first 
marginalize our 4-parameter x 2 analysis with respect to 
the remaining three neutrino parameters. The results are 
shown in Fig. 3. One can see that the new data allow us 
to constraint e, the flavor changing parameter, while for 
the flavor conserving case there is still room for relatively 
large values of e' that correspond to the solution in the 
dark side of the neutrino oscillation parameters. These 
are the bounds we obtained at the 90% C.L. 

-0.41 <e< 0.06 (8) 
-0.50 < e' < 0.19 & 0.89 < e' < 0.99 (9) 

The above limits are in good agreement with the fore- 
cast made in Fig. 3 of [33], assuming the best possible 
determination of the neutrino mixing parameters due to 
KamLAND . In fact, they are even a bit better than ex- 
pected from the improvement of KamLAND data only. 
The reason for this is the subsequent improvement of so- 
lar data, which slightly improved their sensitivity to the 
nonstandard interactions. 

For the flavour-changing effective coupling e, one can 
use the first expression in Eq. (4) to translate the bound 
obtained in Eq. (8) into a limit over the vectorial coupling 
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Figure 3: Bounds on e and e' from the combined analysis of Solar plus KamLAND neutrino data. One sees that the e parameter 
is now constrained to only one region while for the case of e' there are two possible regions, one corresponding to the standard 
light side solution and the other one (less favored) to the so called dark side region. 



-0.08 < ei; < 0.58 (90% C.L.) (10) 

where we have used the best fit value for the atmospheric 
mixing angle [56]. So far, the strongest limit on this 
parameter is \e%*\ < 0.5 [53]. From our analysis, we 
see that solar neutrino data are not only sensitive to the 
sign of this coupling but also we have improved the lower 
bound. 

C. Effects in neutrino detection 

The presence of nonstandard interactions can also af- 
fect the detection processes at some experiments. In par- 
ticular, the cross section for the neutral current detection 
reaction at SNO: 

v + d^ v+p + n (11) 

is proportional to g\. where gA is the coupling of the neu- 
trino current to the axial isovector hadronic current [62]. 
Therefore, the presence of an axial nonstandard coupling 
would give rise to an extra contribution to the NC signal 
at the SNO experiment. This nonstandard contribution 
can be parametrised in the following way [53] : 

4>NC~fB0. + 2e A ), (12) 

where terms of order e\ have been neglected. Here fs 
denotes the boron neutrino flux, and the effective axial 



coupling ea is defined as in [53] : 

e A = - J2 ( P ™) N c4i, (13) 

once the nonstandard axial couplings with up-type 
quarks are set to zero. Note that e^a = £ ta ~ 4m de- 
noting the couplings entering in the effective Lagrangian 
shown in Eq. (1). Thus, e A is independent of the effec- 
tive couplings e and e' defined in Eq. (4) . So far we have 
assumed in our analysis that ea = 0. This assumption is 
well justified due to the good agreement between the SNO 
NC measurement: 4> S N N C ° = 5.54±^;|f (stat) ±%f A (syst) x 
10 6 cm _2 s _1 [11] and the SSM prediction for the boron 
flux f B = 5.94 ± 0.65 x 10 6 cm^s" 1 [48]. However, we 
now relax this assumption by including the effect of the 
new parameter ea- 

The results obtained in a generalized 5-parameter- 
analysis which takes into account the presence of a non- 
zero axial component of the NSI are summarized in Fig. 
4. In the left panel we compare the allowed regions at 
90, 95 and 99% C.L. obtained in the full 5-parameter- 
analysis (filled/colored regions) with the ones obtained 
in the previous section, neglecting the effect of the axial 
NSI couplings (hollow lines). One sees that both analy- 
sis are consistent, though, as expected, the inclusion of 
the axial parameter in the analysis somewhat extends 
the allowed region. In the right panel we show the A\ 2 
profile as a function of the effective axial coupling ea- 
There, one sees that the neutrino data clearly prefers 
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Figure 4: The shaded regions in the left panel show the updated analysis of all solar neutrino data combined with recent 
KamLAND results in the presence of the axial NSI coupling. The lines show, for comparison, the updated regions obtained 
only with the vector NSI couplings. We show in the right panel the x 2 profile with respect to the axial NSI parameter obtained 
using all solar neutrino data combined with KamLAND . One sees how the data prefer ea ~ 0. 



£a ~ 0, thanks to the good agreement between the pre- 
dicted boron neutrino flux and the NC observations at 
SNO, as stated above. We obtain we following allowed 
range at 90% C.L. 

-0.14 < e A < 0.06 (14) 

Using Eq. (13), the above bound on the effective axial 
coupling ea can be translated into individual bounds on 
the NSI parameters ejj~ . Since we are neglecting the non- 
standard interactions of the muon neutrino, this formula 
depends only on the probabilities (P ee ) NC , (P eT ) NC and 
on the NSI couplings e^f and e^ 4 . From the recent values 
for the average probabilities reported by SNO [11]: 

(P ee ) jvc = 0.30 ±0.03 
(Per) NC =0.35 ±0.02 (15) 

one gets an allowed region in the parameter space (e^f, 
£er)i represented (at the 90% C.L.) as a diagonal band 
in Fig. 6 in the next section. There, it will used com- 
bined with neutrino laboratory data to obtain improved 
constraints on the neutrino NSI couplings. 



III. CONSTRAINTS ON NON-UNIVERSAL NSI 

In addition of solar + KamLAND, laboratory experi- 
ments measuring neutrino-nucleon scattering show sen- 
sitivity to neutrino non-standard interactions on d-type- 
quarks. In particular, here we will combine the results 
of the accelerator experiment CHARM together with the 
ones in Sec. II in order to obtain stronger constraints 
on the NSI parameters. Given the sensitivity of the con- 
sidered experiments to different NSI parameters, we have 
been forced to simplify the analyses reducing the number 
of parameters by focusing on the flavour-conserving non- 
universal nonstandard couplings. Within this approxi- 
mation, the parameters relevant for each experiment are 
shown in Table I. 



Data 
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e dv 
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Solar propagation 


/ 


/ 






Solar NC detection 






/ 
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KamLAND propagation 


/ 


/ 






CHARM detection 


/ 




/ 





Table I: Sensitivity of neutrino experiments to flavor conserv- 
ing NSI parameters 
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A. CHARM 

We now turn to the analysis of CHARM data. 
CHARM was an accelerator experiment measuring the 
ratio of the neutral current to the charge current cross 
section for electron (anti)neutrinos off quarks. We have 
used the results reported by the CHARM experiment for 
the v e q — > vq cross section. In particular the experiment 



measured the relation [63] , 

a(u e N -> v e X) + a(D e N —> vX) 



a{v e N -> eX) + a{D e N -> eX) 
{~9Le) 2 + {~9Ref = 0.406 ±0.140 



(16) 



The most general expression for (gL,Re) 2 including all 
types of NSI parameters is given by 



(g Le ) 2 = {gt + etf + E ? + (at + ef e ) 2 + £ l4el (17) 

(9Hef = (ff5 + 4f ) 2 + E l£f I' + + 4f ) 2 + E (18) 

i 



Here we are interested in the flavor conserving d-type 
quark interaction. Then, we will neglect all flavor- 
changing nonstandard contributions, implying that 
e «e = 0, as well as nonstandard couplings to the u-type 
quark, so that e^f = 0. In this case our simplified ex- 
pression for Eqs. (17) and (18) would be 

(hef = (glf + (st + \(ef e + e d J)) 2 (19) 
fee) 2 = (g u R ) 2 + (g d R + \(ef e - e d J)) 2 (20) 

Then, we can compute the x 2 f° r t ne CHARM data: 

( ft e _ r>theo( p dV p dA\\ 2 
2 _ V-"- _ V fc ee i °ee )) /„-. n 

(°rF 

where R e and crfj are defined by the result given in 
Eq. (16). 

The constraints in the (e^J , e d f) plane at 68, 90, 95 
and 99% C.L. obtained from this analysis are shown in 
Fig. 5. 

B. Combined analysis 

In the previous sections we have discussed the sensi- 
tivity of solar experiments, KamLAND and CHARM to 
the nonstandard interactions separately. In this subsec- 
tion we exploit the complementarity of the information 
we can get from the different experiments by combining 
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Figure 5: Constraints at 68, 90, 95 and 99% C.L. on the 
neutrino NSI couplings ef£ and e^t from CHARM data. 

CHARM data with our previous results from the analy- 
sis of the solar and KamLAND data. This enables us to 
obtain stronger constraints on the NSI couplings. 

The regions for the vector (left) and axial-vector (right) 
NSI couplings allowed by the global analysis are given in 
Fig. 6, where they are also compared with the constraints 
coming only from the CHARM data and that from the 
solar plus KamLAND data. First, the bounds obtained 
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Figure 6: Constraints on the vector (left) and axial-vector (right) NSI couplings from our global analysis at 90, 95 and 99% 
C.L., and from the separate solar plus KamLAND and CHARM data sets (dashed lines). 



Table II: New constraints on the vectorial and axial NSI couplings at 90% C.L. obtained from the analysis of CHARM data 
alone and from our global analysis combining CHARM with solar and KamLAND results. 



vectorial couplings 



global -0.5 < e d l < 1.2 -1.8 < e% v T < 4.4 



one parameter at a time 



CHARM -0.5 < eiY < 1.2 

global -0.2 < e A Z < 0.5 -1.1 < e d X < 0.4 & 1.6 < e d X < 2.2 



axial couplings 



global 


-0.4 < eti 


< 1.4 -1.5 < e d 4 


< 0.7 


one parameter at a time 


CHARM 


-0.4 < e d J 


< 1.4 




global 


-0.2 < e d J 


< 0.3 -0.2 < e d 4 


< 0.4 



from the analysis of CHARM data and shown in Fig. 5, 
have been translated into two independent bounds on 
ef£ and e^f (see vertical bands in Fig. 6). The regions 
allowed by the solar+KamLAND combination (diagonal 
bands in the figure) have been derived from the bounds 
on the non-universal nonstandard neutrino interactions 
obtained in Sec. II. In particular, the limits on the vecto- 
rial couplings e^Y an d sfr come from the allowed values 
for the effective coupling e' in Eq. (8), after using the def- 
inition of e' of Eq. (4) and the la allowed region for the 



atmospheric mixing angle 023 [56]. Note the existence 
of two allowed islands, in correspondence with the two 
allowed regions of neutrino oscillation parameters in the 
presence of NSI (see upper-right panel at Fig. 1). The 
lower one corresponds to the usual LMA solution, while 
the upper island comes from the solution in the dark 
side. On the other hand, using the average probabilities 
in Eq. (15) we have reanalysed the results obtained for 
the effective axial coupling ea in Eq. (14). This results 
in a constraint for the axial couplings e d ^ and ef^. From 
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the two panels of Fig. 6 one sees that, as expected, there 
is a degeneracy in the determination of the two vectorial 
and axial parameters from solar and KamLAND data 
only. After the combination with CHARM we break this 
degeneracy and obtain the allowed regions shown in color 
at Fig. 6. 

In Tabic II we quote the 90% C.L. allowed intervals 
for the couplings e^Y, 4n £ tt an d e^r arising from 
the combined analysis, as taken directly from Fig. 6. 
In order to compare with previous bounds obtained in 
a one-parameter-at-a-time analysis [53], we also present 
the results obtained following that approach in our anal- 
ysis of CHARM data alone and in the global analysis in- 
cluding solar and KamLAND data as well. In this case, 
we see how the combination with solar and KamLAND 
data improves significantly the existing bounds on the 
electron- neutrino NSI couplings, obtained using CHARM 
data only [53]. Now, concerning the tau-neutrino NSI 
couplings, we have improved the existing bounds derived 
from the invisible decay width measurement of the Z from 
LEP data: \e*Y\ < 2.7, \e^\ < 1.8 [53]. 

IV. CONCLUSIONS 

We have updated the solar neutrino analysis for the 
case of NSI of neutrinos with d-type quark by includ- 
ing new solar data from SNO phase III and Borexino, as 
well as new KamLAND data and updated solar fluxes. 
We have found that the additional dark-side of neutrino 
parameter space found in Ref. [33] still survives, while 
the previous LMA-0 and LMA-II which were present be- 
fore have now disappeared as a result of the new data. 
The issue arises of how to lift this degeneracy in future 
studies. First we note that, since KamLAND is basically 
insensitive to matter effects, it will not help in resolving 
the degeneracy, as explicitly verified in Fig. 3 of Ref. [33] . 
Next comes improved solar neutrino data. The form of 
the expected neutrino survival probability shown in Fig. 2 



suggests that the best region to discriminate the degen- 
erate solution from the normal one is the intermediate 
energy solar neutrino region, where the relevant experi- 
ments are Borexino [12-14] and KamLAND-solar [64], as 
well as the low energy threshold analysis expected from 
Supcr-K-III and SNO [60, 61]. 

Further information relevant to lift the degeneracy may 
come from atmospheric and laboratory data. Indeed the 
LMA-D solution induced by non-standard interactions 
of neutrinos with quarks may become inconsistent with 
atmospheric and laboratory data [33, 65]. As noted in 
Ref. [33] currently it is not. Should the situation change 
with improved data we note that still the degeneracy will 
not disappear since the NSI couplings may affect not only 
down- type quark species but also up- type quarks and/or 
electrons. The analysis in this case would introduce new 
parameters. Therefore, we conclude that the neutrino 
oscillation interpretation of solar neutrino data is still 
fragile with respect to the presence of nonstandard inter- 
actions. 

In summary, we have studied the limits on the non- 
standard interaction couplings e^Y , ef^, ef£, e d e ^ and 
ef^ from present neutrino data. Thanks to the combina- 
tion of solar and KamLAD neutrino data with the results 
from the CHARM experiment, we have given improved 
bounds on the vector and axial NSI couplings involv- 
ing electron and tau neutrino interactions on down-type 
quarks. 
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